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■ of the Rothwarf- Taylor (RT) equation,^' ^'^^ which is first presented in ref. 10. This type of 
phenomenological equation is claimed to be theoretically derived in refs. 11 and 12 with the 
use of the kinetic equation for the electron distribution function. 

The temporal evolution of the superfluid weight is observed in the experiment with the use 
of an optical pump and a THz probe, and it is found that the photoinduced change in the 
superfluid weight decays nonexponentially. The RT equation has been considered to be suitable 
for describing nonexponential relaxation as the bimolecular relaxation, because this equation 
includes the quadratic term of the nonequilibrium quasiparticle density as the recombination 
of the quasiparticles. However, the quadratic term is not allowed in the perturbation expansion 
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00 , In the transient state of d-wave superconductors, we investigate the temporal variation of 

photoinduced changes in the superfluid weight. We derive the formula that relates the non- 
pH , linear response function to the nonequilibrium distribution function. The latter quantity is 

O , obtained by solving the kinetic equation with the electron-electron and the electron-phonon 

, interaction included. By numerical calculations, a nonexponential decay is found at low tem- 

peratures in contrast to the usual exponential decay at high temperatures. The nonexpo- 
C/5 , nential decay originates from the nonmonotonous temporal variation of the nonequilibrium 

"j^ , distribution function at low energies. The main physical process that causes this behavior is 

^ , not the recombination of quasiparticles as previous phenomenological studies suggested, but 

, the absorption of phonons. 

■ 

Q ' KEYWORDS: nonequilibrium superconductivity, transient state, pump-probe, superfluid 

1 weight, vertex correction, unconventional superconductors 

> ■ 

OO ' 

\C) !• Introduction 
OO 

There have been many optical spectroscopic experiments on nonequilibrium superconduc- 
^> tors recently, and they have usually been performed by pump-probe techniques. ^^^^ As for 
\ the interpretation of experimental results, the photoinduced change in reflectivity is often as- 
sumed to be directly proportional to the excited nonequilibrium electron density. ^''^'^) Then 
the temporal evolution of the photoinduced quasiparticle density is analyzed with the use 
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of the external field for the reason that it does not satisfy the consistency of the equation with 
respect to the order of the pump intensity. Therefore, a question arises about the origin of 
the nonexponential relaxation in the case that the intensity of the pump beam is low. 

In this study with regard to this problem, we reexamine two points that have been usually 
assumed in previous works. One is about the relation between the photoinduccd change in 
optical conductivity and the nonequilibrium quasiparticle density. This relation is not self- 
evident, but it has not been investigated so far, in addition to the expression of the former 
quantity itself. The other is about the temporal evolution of the nonequilibrium quasipar- 
ticle density. The phenomenological description given by the RT equation^*^' is based on 
averaged quasiparticle density. This is also open to doubt for reason that the nonequilibrium 
quasiparticle density should have a dependence on energy. The kinetic equation is not simply 
averaged over energy because of the existence of the interaction effect as the kinetic equation 
in the normal state^^"^^^ indicates. 

To investigate these problems, we microscopically calculate the response function of the 
transient state. Our calculation shows that the vertex correction term is predominant in pho- 
toinduced changes in the superfluid weight, which is same as that in the case of the steady 
state in a previous study. -'^'^^ The nonequilibrium distribution function is related to the re- 
sponse function only through the interaction effect (the vertex correction term). The kinetic 
equation for this function is solved with the use of the electron-electron and electron-phonon 
interactions as the collision integral. The numerical calculation shows that, at low tempera- 
tures, the photoinduccd change in the transient reflectivity is not directly proportional to the 
number of photocxcitcd quasiparticlcs, and the nonexponential decay docs not originate from 
the recombination term, but from the absorption of phonons that enhances the nonequilib- 
rium electron distribution. The effect of nonequilibrium phonons is also considered, and this 
makes the relaxation dynamics slow as a result of the interaction effect. 

2. Response Function in the Transient States 

The current under the pump (AP") and probe (A^^^) beam is written as 

Here, A^^ = A'^^[e-^'^"^^-^°)^ + e-AT2(Q+Qo)V4] and AY^ = e''^U'p,e-^*''^V4 in the case of 
the Gaussian pulse; Af = A^^ cos(J^or)e-(^/^^)' and A?' = A^, cos(wor)e-(^/^*)' {A'^^ = 
and A'^, = A^.^At). Then A^ = Jdr^^^A^ and AS^ = Jdre'^^^A^^ in 
which t is the time delay of the probe pulse after the pump excitation. We put = to 
probe the the superfluid weight. The expression of 

is given in ref. 18, but the rewritten 
form is presented below, which clarifles its relation to the nonequilibrium distribution function 
and its temporal evolution. The integration of K^^\—uj, —Q + uj/2, Q + uj/2) by u and $7 will 
lead to the temporal variation of the superfluid weight, if possible. However it is difficult to 
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perform this integration numerically for the reason that the calculation at a very small energy 
scale is required to determine the variation for a long time. Then we transform the expression 

of 

into a feasible one, as shown below. 

We make the following assumptions to calculate the nonlinear response function. The 
order parameter takes a constant value and is not affected by the external field. The origin of 
the pairing interaction is not specified (the existence of the superconducting gap is assumed 
at the outset). The momentum dependence of the vertex correction is weak, which leads to 
the omission of and S^^) is the self-energy that includes the external fields of the 

n-th order). We take the local approximation as discussed in ref. 17. 

Then the nonlinear response function is written as K^^^ = Nq Jpg J dev^gi^} . Here, is 
the quasiparticle velocity, the summation is transformed into the integration A'^o fpg f d^ 
(A'^o is the density of states at the Fermi surface), and gi"^} = tanh ^{g^^"^^ — gf}'^'^)+ 5cfc '"^^ 
g = f (the integration of Green's function G by the electron dispersion ^), and R, A, and 
(a) indicate the retarded, advanced, and anomalous parts of Green's function, respectively, as 
introduced in ref. 20, in which the microscopic formulation for the nonequilibrium supercon- 
ductors using Green's function is given. (3) indicates the third-order of the external fields, and 
we represent g (g) as Green's function, which includes (does not include) the external fields H^^. 
The predominant term gi^F^ is determined by ^c^fe'^^^ =1^1 ^oif^^ H_n+uj/2Hn+c,/2H-^, 
which satisfies the following equation with reference to Eliasliberg's formulation.^*') 



f dw R(2) -.(a)(1) _ (a)(1) ^A(2) yia){2) A(l) _ R(l) ^(a)(2), 
/ 2^'-^t^~'^9e—u>,e ye,e—u! e—u),e' e,e—(jjye—uj,e "e,€— a; e— a;,eJ" 

Here, = tanh 27^, e^^^^ = e - and H^^ = VkA^^ with the external field A^^. The 

functions g in this equation satisfy similar equations. By taking only the predominant terms 
into account, gi'^]^^^ is written as 



A{2} . 
e—ui,ei 



1 /■dL^r^(a)(2) A(l) _ R(l) y. (a)(2) i „ 



(1) 



Here, 

' ~ eR-{e-n + uj/2)^ eR-{e-n + oj/2)R^eR-{€ + n + uj/2)^ e^-{e + n + oj/2)R' 
To derive the temporal variation of K^^\ we consider the following integration with use of 
some function /^^(ri); 

fdu f di^ ^ ^pX,exp[-tV(At^ + ArV2)] 

(2) 



3/18 



J. Phys. Soc. Jpn. Full Paper 

If the pulse widths At and AT are several femtoseconds, the effective range of frequency is 
of the order of the superconducting gap (Aq). Here, we consider the case that QqAT ^ 1, 
r^o ^ ^0) Ei'iid the dependences of ftoi^) on lo and Q are weak within \uj\ < 1/At and |il=F^o| < 
1/ AT, respectively. Firstly, we examine the temporal variation of the first term of eq. (1). If 
we take account of the above consideration and the weak dependence of this term on oj, we 
can apply the integration eq. (2) to this term. This term is negligible in comparison with the 
vertex correction term discussed below because of the smallness of Ro{Qo), which is the same 
as that in the case of ref. 17, in addition to the exponential factor. Next, we consider the second 
term of eq. (1) (the vertex correction term). By using the relations ^^g+i^ — H-^^dg^^/de and 
dg^^/de = —{dgj}^/de)*, the vertex correction term is written as 

1 f d^r^(a)(2) A(l) m) y{a){2hrr ^. f -n J ^^tA ^{a)(.2) . . rr 

Here, we put ^^^^^2 e-ui/2 ~ ^iTe^^^C'^)' which is a functional of g^l^'^\u), and its functional 
form is determined by specifying the interaction and the self-energy. g^l^'^\i^) satisfies the 
following equation: 

lie + u/2f -{e- u/2)^]gif'^ (c) = J m'Ll^^/2^'nl./2 " (^^^^ " 9t-./2)4T (-) 

~ (.fk,e+uj/2 ~ fk,e-u/2)'^t^}'^\^) + i'^k,e+i^/2 ~ ^k,e-uj/2)fk'}^'^\'^)- 

(3) 

Here, T represents the anomalous (off-diagonal) part of the self-energy, g^^ = —2e^/Z^^ and 



= -2Afe/Z« Z« = sgn(e) J(e«)2 - A^ , and Ru^^, Q) = Rk,e{co, Q) + RkA^, -Q) 



„ , iTe+co/2-Te+n){g^^^^/2~ 9t,e+a) , {Te+n-T^-oj/2){9k^^+^/2- 9k,e+a 
Rk,e[^,^) = TTKr / I r,\A ^ 



{e + oj/2)R-{e + n)A {e-u;/2)R-{e + n)R 

{Te+oj/2 - Te+n){9i^^+n ~ 9t,e-oj/2) C^e+n - T^-uj/2){9k,e+n ~ 9k,e-oj/2) 



{e + n)A- {€-uj/2y {€ + n)R- {e-iu/2y 

This latter quantity is related to the pump-induced term and the initial nonequilibrium dis- 
tribution function, which directly reflects the values of the self-energy. This induced term 
vanishes in the case that the self-energy is absent, as indicated in ref. 17. 

As shown in eq. (3) the dependence of g^^J'^\uj) (and ^[."^''^•'(tj)) on uj is strong, and 
this makes the t-dependence of the vertex correction term different from that of the first 
term of eq. (1) and larger than that. By using the integration f (cij)e^^'^*e~^*^'^^/^ = 

■^7=^/ dr^^"]'^-* (r)e^(*~'^^^/^*^ ~ g^kl^'^\'^)^ which is applicable to the time range t » At, 
we can transform eq. (3) into the following kinetic equation, which describes the long time 



4/18 



J. Phys. Soc. Jpn. Full Paper 

behavior of g'^^}'^^ it): 

Here, a(t) = "^pu^Pt * /(^* +^'^ _ obtain K^^\ the solution of this kinetic equation is 
substituted into in the following vertex correction term: 



FS 



The kinetic equation for nonequilibrium phonons is similarly written as 



UJ 



dt 



(6) 



Here, D^'^^^ and H,^'^"^^ are the retarded (advanced) phonon Green's function and the self- 
energy by phonon-electron interaction, respectively, (a) and (2) indicate the anomalous part 
and the order of the external field, respectively, which is same as that in the case of electrons. 

3. Kinetic Equation of the Nonequilibrium Distribution Function 

In this section, we derive the kinetic equations for the distribution functions of electrons 
and phonons. We put the deviation of the distribution function from the equilibrium state 
by the pump excitation as 5ne{t) and SN^^{t) for electrons and phonons, respectively. We 
consider a two-dimensional system and replace Jpg by J Vk = vp cos ip (vp is the Fermi 
velocity) and = Aq cos 2ip for d-wave superconductors. 

We adopt the second-order perturbation expansion for the electron-electron interaction 
and the one-loop approximation for the electron-phonon interaction as the self-energy. We 
rewrite the nonequilibrium Green function as g^p}J"^\t) = —2Sne{t){g^^^ — g^^^), f!p^J^'^\t) = 
-2SnS){f^,e - f^,e\ and Dl^^^''\t) = 25N^{t){D^ - D^) with the use of - = 
—'it\.uj[5{oj — Uff,) + 6{oj + cjfj))] (for example, see ref. 21). We consider acoustic phonons and 
put LJ^ = Vskplcf)] {vs is the sound velocity). Then the kinetic equation for electrons, eq. (4), 
is written as 

d6n,{t) 1 a{t) fdip^-.^^. 1 vr , . /"^^ d</) 



1 nU^Ni /•/•/• dydyidy2 ff\ , ~, rel-el rr i 

^2^ JJJ (27r)3 JJ deide2(^^,^„^,4,^;,„^,;,„^JH 



Here, ■= - J ^Reg^{e) and S^,^-,^^^ := S[l - cos{(p - ipi) + cos{(p - (^2) - cos(¥3i - (^2)] 
is a delta function. U is the effective short-range Coulomb repulsion energy, which is called 
the on-site Coulomb interaction in the Hubbard model, g is the coefficient of the electron- 
phonon interaction, which is usually written ss g = ^/h/2MiVsVi (Mi is the mass of the 
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ion and Vi is the rcnormalizcd ion potential). In the next section, we use UNq and g'^N^ as 
dimensionless parameters that characterize the strengths of the electron-electron and electron- 
phonon interactions, respectively. The collision terms for the electron-electron and electron- 
phonon interactions are written as 

„l_el 1 ^ Reg^.eReffyi ,ei ^^9^2,^2 ^^gy-ipi+iP2,e-ei+e2 

6,<^;6i,<^i;62,c^2l J cosh cosh |L cosh cosh ^^f^ 

X [(cosh ^f5ne{t) - (cosh ^fSUe^it) + (cosh ^fSUe^it) - (cosh ^~2T^^ )^'^"-6-6i+62 

and 

I^l~j\5n, 5N] = (coth ^ + tanh ^) 9^Jne{t) + (coth ^ - tanh ^) Q^e^^S) 

- (coth ^ - tanh i§j:) g;^ Jne-oj^{t) - (coth ^ + tanh 27^) 5fJ,(5ne+a;^(t) 
(^tanh ^ - tanh g^JN^{t) + (i&nh ^ - tanh 27.) g^JN^{t). 

Here, := / g [Re<7^_^(e T c^</,)Re5^(e) - Ref^-^e T a;0)Re4^(e)^ . 
The kinetic equation for nonequilibrium phonons, eq. (6), is written as 

Here, we add a phenomenological term, jesc^^^u^ (damping by phonon escape), which describes 
the equilibration between the electron-phonon system and the reservoir. 

if'lSrie, SN^] = (coth ^ + tanh ^) g^Jn,{t) - (coth ^ - tanh ^) gX^Sn^it) 

- (coth ^ - tanh ^) g~ Jne-oj^{t) + (coth ^ -|- tanh ^) g+Jne+^^{t) 
+ (tanh ^ - tanh g- JN^{t) - (tanh - tanh 25.) g+JN^{t). 

The energy conservation is also discussed using the kinetic equations. The additional 
energy induced by the external field in the system of electrons and phonons is given by 

2^EkSnE,{t)+Y,^gSN^^{t) = 2No [ de e6nS)+^ f^" ^27rvFkF\ct)\u^^5N^{t). 



k q Y " "V 

With the use of the kinetic equation, it is shown that this quantity is equal to the energy 
injected by the external field, which is written as 



In the above discussion, we put 7esc = for the energy conservation. In the case of 7esc 7^ 0, 
the energy of the electron-phonon system dissipates into the reservoir system. 

The solution of the kinetic equations is related to the nonlinear response function as 
follows. 6ne{t) and 6N^{t) obtained by solving the kinetic equations is substituted to S^'^g*-^-' {t) 
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of the previous section as = S^^e*^' '^\t) + ^'^(*)- Here, 

^' ~ VFkF J J (2vr)2 J J ^^^^ _^ ^^^^ ^ ^^^^ ^ ^^^^ 

X [(cosh §.)^^nei(i) - (cosh ^fSne^{t) + (cosh '"ar '' )^(^ne-ei+e2(^)] , 

and 

5.(a)ei-ph^^) = - TTiNog'' J'^^ ^c^^Reg^.^,,.,^ [(coth ^ - tanh ^) 5n,.^^{t) - (tanh ^ - tanh 5^) ^iV^t)] 

+ Re5^_^^,+^^[(coth ^ + tanh 27.) 5ne+u^^{t) - (|tanh ^ - tanh 5N^{t)\}. 
Then this S^"e'~^'* (t) is substituted into eq. (5) , which gives the nonhnear response function. 
4. Results 

The results of the numerical calculation are shown below. The superconducting gap Aq is 
taken as the unit of energy Aq = 1.0 (this leads to the superconducting transition temperature 
Tc = 0.465), and we put < Vsk^- We fix the values of several parameters as follows: 

= 6.0,22) ^^^p ^ 2.0372, Vs/v-p = 0.05 (which leads to the Fermi energy Ep = v^k^/^ ~ 20), 
UNq = 0.2, and g^No = 0.05 (this value corresponds to A = 0.1 as the electron-phonon 
coupling constant). We add a small S = —0.01 to the imaginary part of the self-energy as 
effective impurities and a finite mean free path. 

In the calculation of the kinetic equation eq. (7), the following approximation is 
used to make a multiple integration in the electron-electron collision term a feasible one: 

We perform the same approximation for This leads to violation of the momentum 

conservation, although the energy conservation is satisfied. The numerical calculation indicates 
that, by this replacement, the interaction effect is underestimated at low energies (|e| < Aq) 
owing to averaging the angular dependence. However, the quantitative difference is small es- 
pecially at |e| > Aq. Then it is not considered that this approximation causes qualitative 
changes for numerical results. 

Below we consider two typical cases for a nonequilibrium phonon system. The numerical 
calculation of K^^^ shows that calculations with large values of the damping 7esc lead to results 
similar to those in the case of 6N = (phonons in thermal equilibrium). On the other hand, 
the results with small values of 7esc are approximated by those for 7esc = 0. Therefore, we 
show the calculated results only for SN = and 5N ^ with 7esc = 0.0 as characteristic cases 
of the phonon system. 

To solve the kinetic equations, we consider the square pulse and put a(t) = 1 for < t < 
0.25 and a{t) = otherwise. (The long time behavior is not affected if we use a Gaussian 
pulse.) Then we redefine K^^^ as K^^^ = -4c\t)/A'p, from eq. (5) with the use of this a{t), 
and this is the photoinduced change in the superfiuid weight. (We can put J^^\t) 2± Jyc\t) 
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Fig. 1. Dependences of I/K'^^^ on time for various values of T/T^ (numbers in the figure). The 
phonons are supposed to be (a) in thermal equilibrium (5N = 0), and (b) in nonequilibrium state 
{6N ^ 0). 



as noted above.) If we use only S^")^^ (S^")^' p**) as S^")^^) in this equation, we write this 
quantity as itC(3)ei-ei (^(3)ei-ph) ^^^^^ 

The temporal variation of is shown in Fig. 1. If we put Aq = 30 meV, t = 1000 

corresponds to about 22 ps. This is comparable to the range of time taken in the experiments,^^ 
and we present the numerical results for this range hereafter. K^^^ shows a nonexponential 
decay for low T/T^, and it becomes an exponential decay for high T/T^ as in the experimental 
result.^) Here, the nonexponential decay and exponential decay indicate oc 1 + 7* and 

oc e'''*, respectively. At first sight, the former seems to be an approximation of the latter 
using e'"'* 2± 1 + 7* with 7* ^ 1, but this is not the case. For example, at T/T^ = 0.05, 
is fitted to 0.1(1 + 0.0076t) (0.0076t » 1 for large t). Therefore, some qualitative differences 
exist between low and high T/T^. It is shown below that this difference in relaxation dynamics 
arises from the temporal variation of the nonequilibrium quasiparticles ((Jn^), and this traces 
back to the change of the predominant physical process in the collision integral. In the case of 
6N 7^ 0, the relaxation becomes slower than that in the case of 5N = 0. As shown below this 
originates from the existence of 5N ^ which directly affects K^^^ through the electron-phonon 
interaction, in addition to the slow relaxation of nonequilibrium electrons. If we vary g'^No, 
we obtain qualitatively the same results with regard to the t- and T-dependences of 1/K^^\ 
but there are several quantitative differences. In the case of SN = 1/ K^^^ increases, and the 
exponential decay of is observed at lower T/T^ with increasing ^^A^o- This is because 

the equilibration between electrons and phonons becomes faster, and Sn rapidly decreases. 
On the other hand, in the case of SN 7^ 0, the absolute value of decreases with 

increasing g^N^. In this case, the nonequilibrium electrons remain finite, and the interaction 
effect directly enhances K^^\ 
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Fig. 2. Dependences of ii:(3)oi-ph on time for T/T^ ^ 0.1. 5N ^ Q miA 5N ^ indicate the results 
of iir(3)°'~P'^ in the cases of 5N = and 5N ^ 0, respectively. [5n\ and [6N] are components of 
/;i'(3)ei-ph £qj, ^ and include only 5n and 6N in I]^"e°' , respectively. 




400 800 1200 1600 2000 400 800 1200 1600 2000 

t t 

Fig. 3. Dependences of {5n) on time for various values of T/Tc (numbers in the figure). The phonons 
are supposed to be (a) in thermal equilibrium {5N = 0) and (b) in a nonequilibrium state {5N ^ 0). 

To clarify the effect of nonequilibrium phonons on K^'^^ , we decompose iif (^)^'~p'' to con- 
tributions from nonequilibrium electrons and phonons by taking only the Sn term or 6N term 
in T,^^e^~'^^. The temporal variations of these terms are shown in Fig. 2. This result indicates 
that the 6N term has a sufficient contribution to K^^") at a large t. This makes the relaxation 
slower than that in the case of SN = 0. 

Here, we examine the question as to whether the nonlinear response is proportional to 
the nonequilibrium electron density. The temporal variation of the integrated nonequilibrium 
electron density, which is written as {Sne{t)) = de Jpg —===5nf:{t), is shown in Fig. 3. 
In the case of SN = 0, {Sn) shows the exponential decay for high T/Tc, as K^^^ does. For 
low T/Tc, (Sn) does not show the same t-dependence as K^^\ Therefore, in the region of 
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Fig. 4. Dependences of X*^'^'"' °' (ee) and if^'''"' (ep) on time for various values of T/T^. The 
phonons are supposed to be (a) in thermal equilibrium {5N = 0) and (b) in a nonequilibrium 
state {5N ^ 0). 

the nonexponential decay, there is no proportionality between the nonlinear response and 
nonequilibrium electron density, which is different from the assumption adopted in previous 
studies. The similarity in relaxation dynamics between these two quantities is limited to the 
region of the exponential decay. In the case of 6N ^ 0, {6n) does not decrease to 0, but it 
varies toward the nonequilibrium steady state. This state is rapidly achieved for high T/T^. as 
for {6n). The increase of {6n) with t at low T/T^ is understood by considering the temporal 
evolution of 5n^{t) (shown in Fig. 7) and the collision integral I^{t) (in Figs. 6 and 9). As 
discussed there, the distribution of nonequilibrium electrons is enhanced at low energies by 
the absorption of phonons. (The particle conservation is satisfied because (5n_£ = —5ne in our 
calculation.) 

Next, we investigate which of the electron-electron and electron-phonon interactions pre- 
dominate in the nonlinear response. The temporal variations of K^^^^^~'^^ and if^^)^'"^' are 
shown in Fig. 4. For small (large) t, K^^^^^~^^ (7if(3)ei-ph'j -g predominant in K^'^\ This change of 
the predominant term originates from the temporal variation of the functional form of Sne{t). 
Sne{t) has a broad spectrum as a function of e at small t, and its spectrum concentrates at 
low energies as time passes, as discussed below. The time at which K^^'>^^~^^ and K^^^^^^p^ 
intersect is weakly dependent on whether phonons are in the equilibrium or nonequilibrium 
state, which indicates that the effect of nonequilibrium phonons is small in this time range. 
Although K^^'>^^~'P^ decreases with increasing T/T^, K^^'>^^~^^ shows a different behavior at 
large t, which also reflects the e-dependence of Sne{t), but this behavior is not important 
because K^^>^-^^ < i,f(3)ei-ph ^^^^ ^-^^ ^^^^e. 

Hereafter, we consider the microscopic quantities that cause the above behavior of the 
nonlinear response. Firstly, we show the temporal variation of the collision integral, which is 
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Fig. 5. Collision integrals I^{t) for various values of t. T/Tc = 0.1. The phonons are supposed to be 
in thermal equilibrium {dN = 0). The results at T/Tc = 0.3 are shown in the inset. 



the second and third terms on the right side of eq. (7) and written as Ie(t) = lf~^^(t)+if~^^(t) 
hereafter. As indicated above in the cases of K^^^^^"^^ and k(3)^i~p^, in the coUision integral 
Ie{t), the electron-electron interaction predominates the relaxation dynamics over the electron- 
phonon interaction (roughly Ig{t) 2± lf~^^{t)) at small t, and it is the other way around 
{Ie{t) 2± lf~^^(t)) at large t. This holds irrespective of the values of T/Tc and the existence of 
SN. We show the collision integral Ie{t) (only for e > because of I-e{t) = —Ie{t)) at small 
t in Fig. 5 and at large t in Fig. 6. For small t, there seems to be no qualitative difference in 
the e-dependence of /e between T/Tc = 0-1 and 0.3. The energy range in which takes finite 
values is broad, and it takes negative and positive values at high and low energies, respectively. 
This brings about a shift in the weight of the nonequilibrium distribution function from high 
energy to low energy. The results in the case of SN ^ are omitted here, but they are similar 
to those in the case of SN = 0. In contrast to that in the case of small t, e-dependences of 1^ 
are different depending on the values of T/Tc and SN for large t. For T/Tc = 0.3, 1^ becomes 
negative all over the range of e > 0. This property is the same in the case of SN ^ 0, although 
its degree is small. For T/Tc = 0.1, there remains a positive part in /g at low energy. The 
shift in the weight of the nonequilibrium distribution function from high energy to low energy 
occurs, as in the case of small t, but its energy scale becomes narrower. The result for SN ^ 
shows that the negative part in is small and that the decreasing rate of the positive part is 
slower than that of SN = 0. 

The nonequilibrium distribution functions for electrons, Sne{t), at various times t for 
T/Tc = 0-1 and 0.3 in the case of SN = are shown in Fig. 7. (We consider the particle-hole 
symmetric case: Sri-^ = —Srie.) At small t, immediately after the pump excitation, Sne{t) has 
a broad spectrum. The spectrum at high energy rapidly decreases with time. This is because 
the damping effect is large at high energy owing to the electron-electron interaction. Then 
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Fig. 6. Collision integrals I,{t) for various values of t. (a) T/T^ = 0.1 and (b) T/Tc = 0.3. The 
phonons are supposed to be in thermal equilibrium {5N = 0). The results in the case of 6N ^ 
are shown in the inset. 
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Fig. 7. Nonequilibrium distribution function 6ne{t) (a) at small values of t and (b) at large values of 
t. T/Tc = 0.1. The results for T/Tc = 0.3 are shown in the inset. The phonons are supposed to be 
in thermal equilibrium {6N = 0). 



SnJt) takes a functional form, which is similar to that of the r*-model:^^) .J* , .1 , ~ 

[2Tcosh(7^r)]'^ ' where T* is the temperature that characterizes nonequilibrium electrons and 
T* > T. For T/Tc = 0.1, the shift of the spectrum from high energy to low energy occurs. 
(This indicates a nonthermal state, as discussed in ref. 14, with taking account of a change in 
the sign of d5n^/dt ~ j^-'p^ fQj- large t.) On the other hand, 5n^{t) starts to decrease all over 
the range of e > after a certain t for T/Tc = 0.3. This behavior is understood by examining 
the e-dependences of the collision integral, as shown above. For 5N ^ 0, these tendencies are 
qualitatively the same, which is also indicated in the results of the collision integral. 

The energy injected by a pump beam is transferred to the phonon system via the electron 
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Fig. 8. io6Nu{t) (a) at small values of t and (b) at large values of t. T/T^ = 0.1. 



system. In the case of 5N = 0, we presume that this energy transferred to the phonon system 
dissipates by the damping effect 7esc and that there is no influence on the electron system 
by SN. On the other hand, SN is finite if the effect of 7esc is small. The distributions of 
nonequilibrium phonons u!SN^{t) at various values of t for T/Tc = 0.1 are shown in Fig. 8. 
The tendency in the changes of the spectrum is similar to that observed in nonequilibrium 
electrons. However, the decrease in the spectrum does not occur in a phonon system, and 
nonequilibrium phonons accumulate at small cj with increasing t. For high T/Tc, SN shows a 
similar but broader spectrum than that for low T/Tc. 

Finally, we investigate which of the physical processes causes the nonexponential decay of 
K^^^ by increasing Sn^. We rewrite the collision integral by the electron-phonon interaction 
jei-ph ._ _j^n]\j^^g2 J^^^ ^Qj^j^~P^^Sne, SN^] to Specify the physical processes in this term. 
In the case of e > 0, is decomposed into three terms, ~ -^^ e + -^^ e + -^^ e' which 

each term is written as follows. 

lie = 25^,e'^^'^K(l - ^e-c.,)(l + N^,) " u^.^^il - n,)iV,J^(e - cv^). 
lie = -2</^'^[ne+a,,(l - n,)(l + N^,) - n,(l - n,+,,)iV,J. 

4,e = '^ale^^^^nen^^-eil + N^^) - (1 - n^^.,){l - n,)N^^]e{u^ - e). 

Here, (J^^) operates n or N ((2) indicates the order of the external fields), and then 
Si^^Ue = Sne{t) and S^^'^N^^^ = SN^^{t) (other n and N are replaced by n° = l/(e^/^ + 1) 
and N^^ = l/^e"*/'^ — 1), respectively). Each term describes the emission of phonons (/"), 
the absorption of phonons (J^), and the recombination term (J'^), respectively. (The recom- 
bination does not mean that quasiparticles recombine into Cooper pairs, as the misleading 
picture in ref. 11 suggests. This term exists in the case of the electron-phonon interaction for 
d-wave superconductors.) We substitute these terms into the above collision integral lf~^^, 
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Fig. 9. Energy dependences of if-P^''{t) and at t = 1500 for (a) T/T^ = 0.1 and (b) T/T^ = 

0.3. The phonons are supposed to be in thermal equihbrium {5N = 0). 



and write each term as 7"'^''^. 

The energy dependences of if~^^{t) and 7^'^''^(t) in the case of 6N = are shown in Fig. 9. 
A scattering term (emission of phonons) /" and the recombination term are negative, and 
another scattering term (absorption of phonons) is positive. At low T/Tc, l\ is predominant 
over 7" + 7^ , and then the cohision term takes a positive value at low energy. Although 

the absolute values of each term decrease with increasing T/Tc, jf~^^ takes a negative value at 
r/Tc = 0.3 because of a rapid decrease in the degree of absorption of phonons. The negative 
jei-ph fQj. g > Q brings about a decrease in 5n^{t) all over e > 0, as shown above. This 
causes the exponential decay of K^^\ On the other hand, the positive 7|^~p^ at low energy 
increases 5ne{t), and then reduces the decreasing rate of K^^^ as compared to that of the 
exponential decay. Therefore, we regard the absorption of phonons as the main process for the 
nonexponential decay of K^^^ . This is in contrast to the phenomenological interpretation based 
on the RT equation,^) in which the recombination term causes the nonexponential relaxation 
as the bimolecular decay. In the case of 5N ^ 0, the energy dependences of if~'^^{t) and 
ja,6,c^^-j are qualitatively similar to those of 5N = 0. However, in this case, the relaxation 
dynamics is slower than that of SN = because of the presence 
above. 

This predominance of the phonon-absorption term at low temperatures is restricted within 
some range of time. By performing the calculation further at large t, it is found that the 
recombination term becomes predominant over the absorption of phonons. At T/T^ = 0.1, 
for instance, lf~'^^{t) takes negative values in the entire range of e > at approximately 
t 2± 5000. Then, at about this time, deviates from the t-linear behavior, and ceases to 

show the nonexponential decay. 
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In this study, we investigated the photoinduced change in the superfluid weight in the 
transient state of d-wave superconductors. We clarified how the nonexponential decay occurs. 
At first, we obtained the formula that relates the nonequilibrium distribution function to the 
physical response function. Although the former quantity is usually discussed in theoretical 
works, it is the latter quantity that is necessary for a comparison with the experiments. In 
this derivation, we found that the vertex correction is predominant and that the nonequi- 
librium distribution function is effective only through this interaction term. We numerically 
solved the kinetic equation for the nonequilibrium distribution function with taking account 
of the electron-electron and electron-phonon interactions and substituted its solution into the 
expression of the nonlinear response function. 

The numerical calculation shows that the electron-phonon interaction predominates over 
the electron-electron interaction in the long-time behavior. In contrast to the previous studies 
based on the phenomenological RT equation, the nonexponential decay does not originate 
from the bimolecular recombination. Rather, it results from the enhancement of the nonequi- 
librium distribution at low energies, which is caused by the absorption of phonons. This fact 
is revealed using the nonequilibrium distribution functions that couple with each other at 
different energies through the interaction effect, and is not known from the RT equation in 
which the quasiparticle density averaged over energy is used. This leads to an explanation 
of the nonexponential decay which is consistent in terms of the order of the external field 
under the condition that the pumping intensity is low. (As noted in ref. 17, this condition is 
presumably satisfied in the experiment^^ with reference to the excitation fluence in ref. 7.) 

Finally, we comment on several problems related to but beyond the scope of this paper. 
Strictly speaking, it requires numerical integration with a very small energy scale to discuss 
the long-time behavior of a response function. In this paper, however, we avoided this by 
performing Fourier transformation at the outset. This is achieved at the cost of accuracy 
in the short-time behavior of a response function. Therefore, our formulation is not suitable 
for discussing physical quantities such as the exact form of the spectrum at this time scale; 
however, the results presented in this paper are not considered to be affected by these fine 
structures because of the rapid smoothing by the electron-electron interaction. 

There are very few experiments that probed low-energy phenomena such as the superfluid 
weight reported in ref. 5. Most pump-probe experiments have been performed with the use 
of an optical probe beam. The formulation in this paper is restricted to the case in which 
the probe frequency is zero. The extension to a finite probe frequency is required to discuss 
optical probe cases, but it will not alter the importance of the interaction effect if we consider 
the case of the local limit. 

The formulation based on the local limit is appropriate for cuprate superconductors as 
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discussed in ref. 17. Other superconductors are also investigated by nonlinear optical spec- 
troscopy, and the pump-probe experiments have been performed in iron-based superconduc- 
tors (for example, ref. 24.) Infrared spectroscopy indicates that the nonlocal limit is realized 
in superconductors of this kind,^^) and in this case the Mattis-Bardeen formula^^-* is valid. 
It is possible that the nonlocal limit gives a different result regarding the nonlinear response 
from the local limit in which the interaction effect is essential for the optical response. 
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